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Abstract 

We construct a random matrix model that, in the large N limit, reduces to the low 
energy limit of the QCD partition function put forward by Leutwyler and Smilga. This 
equivalence holds for an arbitrary number of flavors and any value of the QCD vacuum 
angle. In this model, moments of the inverse squares of the eigenvalues of the Dirac 
operator obey sum rules, which we conjecture to be universal. In other words, the validity 
of the sum rules depends only on the symmetries of the theory but not on its details. To 
illustrate this point we show that the sum rules hold for an interacting liquid of instantons. 
The physical interpretation is that the way the thermodynamic limit of the spectral density 
near zero is approached is universal. However, its value, i.e. the chiral condensate, is not. 
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1. Introduction 

Since our main understanding of nonperturbative phenomena in QCD comes from 
numerical simulations analytical results in this direction are most welcome [0, @]. In 
particular the issue of quark confinement has not been resolved. What is better known 
is the mechanism of chiral symmetry breaking. In the chiral limit the QCD partition 
function is invariant under SU{Nf) x SU {Nf) x f/(l), whereas the axial U{1) symmetry is 
broken explicitly by quantum fluctuations. For more than one flavor {Nf > 1) the chiral 
symmetry is spontaneously broken down to SU{Nf) x f/(l) by the formation of quark 
condensates. According to Goldstone's theorem this leads to Nj — 1 massless excitations. 
At low energy, these are the relevant degrees of freedom of the QCD partition function, and 
it is possible to write down an effective Lagrangian for an arbitrary value of the vacuum 
angle 9 in terms of these degrees of freedom only [^, |^, |^. The partition function can 
be simplified further by suppressing the space-time dependence of the Goldstone modes 
altogether. 

Recently, Leutwyler and Smilga evaluated the QCD partition function in this limit. 
This enabled them to write down sum rules for the spectrum of the Dirac operator. What 
enters in these sum rules is what we call the microscopic spectral density as opposed to the 
continuum spectral density that enters in the calculation of the condensate p . The latter 
one is related to the thermodynamic limit of the spectral density, whereas the former one 
provides information on the way the thermodynamic limit is approached. The question 
we would like to address is in how far the sum rules are specific for the QCD partition 
function, and in what sense they are universal. Since the partition function of 0| only 
involves constant fields one would suspect that the detailed structure of the QCD vacuum 
is not important, and that the sum rules only depend on the symmetries of the theory. 

This adagio is taken to the extreme by constructing a model with the correct chiral 
structure, but that apart from this does not contain any other information. As is well 
known from random matrix theory, the minimization of information leads to gaussian 
random matrix ensembles [0. The appropriate ensemble will be constructed in section 2, 
and with the help of mathematical techniques developed in the framework of Anderson 
localization and compound nucleus scattering ^ we are able to derive the partition 
function that was obtained by Leutwyler and Smilga (see section 3). Sum rules for one 
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and two flavors are derived in section 4. The structure of the random matrix model 
was inspired by the semiclassical approximation to the QCD partition function where 
instantons are the main degrees of freedom llO, [Tl|, IT^. Therefore, it is natural to check 



the validity of these sum rules for an instanton liquid model ||T3[ of the QCD vacuum, a 
task which is carried out in section 5. Concluding remarks are made in section 6. Finally, 
in appendix A sum rules for a finite number of degrees of freedom are derived. 



2. Random matrix model 

The Euclidean QCD partition function for Nf flavors and nonzero vacuum angle 6 can 
be written as 

Nf 

^QCD=Ee^'^< nn'(An+^;K>5.(^), (2.1) 

u f=l n 

where the product is over the positive eigenvalues of the Dirac operator, and the aver- 
age < ■ ■ ■ >s^{A) is over gauge field configurations with topological quantum number p 
weighted by the gauge field action Sy{A). The topological part of the action, expiOu, 
has been displayed explicitly. Here and below, the mass matrix is taken to be diago- 
nal. According the the Atiyah-Singer theorem, the Dirac equation in a gauge field with 
topological quantum number u has exactly u zero eigenvalues. However, configurations 
with zero total topological charge may actually be composed of spatially well-separated 
components with opposite topological charge. Such configurations will give rise to al- 
most zero modes and thus will play an important role in the chiral dynamics of the QCD 
partition function. Our starting point is that the chiral properties of the QCD partition 
function are determined by zero modes and almost zero modes only. Although it is not 
necessary for the sake of the argument, it is instructive to think of the field configurations 
as a superposition of instantons and N_ anti-instantons. Each isolated instanton or 
anti- instanton has exactly one fermionic zero mode with a definite chirality [|T^. In total 
we have u = A+ — A_ exact zero modes. At finite separation of instantons and anti- 
instantons the remaining modes are no longer exact zero modes of the Dirac equation and 



give rise to nonzero overlap matrix elements T of the Dirac operator |15 



A model describing the zero mode part of the QCD partition function with N zero 



modes or almost zero modes in the Euclidean volume V is defined by 

Z{e) = E /.(iV) E ( e^'^^--^'^ J VTP{T) n det ( ) . (2.2) 

Zero modes of each chirality are treated as independently distributed identical particles, 
hence the binomial factor. The distribution function fi{N) of the total number of zero 
modes = + is peaked at some average value with a width that is much smaller 
than A^. As long as a large number of different values of A^ contribute to the partition 
function with roughly equal probability^, our results do not depend on the detailed shape 
of /i(A^). The average over all gauge field configurations in eq. (2.1) is replaced by an 
average over gaussian distributed overlap matrix elements with distribution function given 

by 

P(T)=exp(-^TTt). (2.3) 

The integration measure VT is the Haar measure. The structure of the overlap matrix, 
with off-diagonal blocks T and and diagonal blocks equal to the quark masses rrif times 
the identity, is dictated by the chirality of the zero modes. The matrix T is a A^+ x A^_ 
matrix, and for zero quark masses, the total overlap matrix has |A^+ — A^_| exact zero 
eigenvalues. The density N/V of the total number of zero modes is kept fixed. As in 
the instanton liquid approximation to the QCD partition function |jl6|], it is considered 



to be an external parameter. The thermodynamic limit can thus be taken by letting 
N ^ oo, where here and below we will omit the bar. A similar model with A^. = A^_ was 



considered in 17, 18 



The order parameter in the study of chiral symmetry breaking is the quark condensate 

< QfQf > defined by 

< QfQf >= lim lim -l-AiogZ(e), (2.4) 

where the order of the limits should be taken as indicated in the formula. In general, 

< qfqf > depends on 6. Its value at = in the limit where for all flavors rrif — > is 
approached from above is denoted by — S. By writing the determinant as the product 



^If ^J.{N) - S{N ~ N) the parity of N+ - is the same as the parity of N = N+ + N^, and the 
partition function has the additional symmetry Z{9 = 0) = Z{9 = n). 



n'(A^ + "^/) one obtains the Banks-Casher formula for S 

S = 7r<pc(0)>z(e=o), (2.5) 

where the average < ■ ■ ■ >z{e) is with respect to the partition function (2.2) (or (2.1) in 
the case of QCD). The continuum spectral density is defined by 

pc{X)= lim, lim lp(A), (2.6) 

all mf 10 W— ►oo iV 

and the spectral density p(A) is 

p{X)=J26{X-X^), (2.7) 

where the eigenvalues ±A„ are the nonzero eigenvalues of the overlap matrix in the chiral 
limit. 

We will also consider a different limit of the derivative with respect of m of the partition 
function, namely, 

1 # 

hm logZ(^). (2.8) 

iv^oo JyP dm J 

Nrrij fixed 

Let us consider the case p = 1 in more detail. Again writing the determinant as a product 
over eigenvalues, one finds 

Nnif fixed Nmf fixed 

(2.9) 

What enters in this expression is what we call the microscopic spectral density defined by 

i!!.". W"^^^' (2.10) 

Nm f fixed 

as opposed to the continuum spectral density defined by in eq. (2.6). In this limit, the 
spectral density function near zero is enlarged proportional to the size (given by A^) of 
the system. Note that ps{x) depends on Nrrif. 



5 



3. Calculation of the partition function 



In order to evaluate the partition function (2.2) the determinant is written as an 
integral over Grassmann variables 



/ 



(3.1) 



where the measure of the Grassmann integration is as usual 

Vijf = Y[d^{dijt, (3.2) 

i 

and the conjugation * is the conjugation of the second kind (i.e., ip** = HI- The 
integral over T is gaussian and can be performed easily. In the partition function this 
results in the factor 

2A^ 

exp— V^/*V^f0f (3.3) 

which represents a 4-fermion interaction. 

The quartic term can be written as a sum of two squares 

v^f>f0f 0f = + 0f>f)(^f v^/ + 0f 0j) - ^(v^/>f - 0fvr)(cv'/ - 0f 

(3.4) 

Each of the two squares can be linearized with the help of a Hubbard-Stratonovich trans- 
formation [0. This allows us to perform the Grassmann integrations at the expense of 
the introduction of the new real valued integration variables a^^ and a-^^ , respectively. 
Apart from an irrelevant overall constant, the partition function reduces to 

Z = V /i(Ar) V f ^ ) e^'(^+-^-) / VaVa det^+ (a + icr + m)det^- {a - icx + m) 



N 

X exp — -Tr (cr + 20") (a — icr). (3.5) 

2A 



As always, the measure of the integral over the matrices a and a is the Haar measure. 
The diagonal mass matrix is denoted by m. 

The complex matrix a + ia can be decomposed in 'polar coordinates' as [TU 

a + ia = UW-\ (3.6) 



where U and V are unitary matrices and A is a diagonal real positive definite matrix. Since 
the r.h.s has Nf more degrees of freedom than the l.h.s., one has to impose constraints 
on the new integration variables. This can be achieved [|19| by restricting U to the coset 
U{Nf)/U{l)^f , where U{1)^^ is the diagonal subgroup of U{Nf). In terms of the new 
variables the partition function reads 

N 



N N. 



X det^'+iUAV-' + me''/''f)det^-{VAU-^ + me~''^''f) exp{~^TiA^), (3.7) 

2a 



where the integral over U is over U{Nf)/U{l)^f and the integral over V is over U{Nf). 
A phase factor exp{—i6Nf) has been absorbed in the integration over V. 

For Nf flavors we have Nf condensates which break down the symmetry of the action 
to U{Nf)/U{l)^f leaving us with Nj-1 Goldstone modes for ^ N^. When the total 
topological charge is zero the phase of the determinant also cancels which provides us 
with an additional zero mode. 

The term proportional to m plays the role of a small symmetry breaking term. The 
integrals over the nonzero modes will be performed by a saddle point integration at 
m = 0, whereas the integrals over the soft modes will be accounted for exactly at a fixed 
value of mN. There are two types of nonzero modes, the phase exp ia of the determinant 
det VU~^ for A^_|_ ^ N_, and the eigenvalues A. The partition function at m = factorizes 
accordingly, 

da{exp{ia) + exp{-ia))^ / J(A)c/Adet^Aexp( r^TrA^). 

(3.8) 

The leading order contribution in 1/A^ of the integral over A can be obtained by a saddle 
point approximation. The saddle point equations for the A integrals read 

Ai = ±A. (3.9) 

The negative solution is not inside the integration manifold and can be omitted. 

The integral over a, which ranges from until 27r , can be executed either before or 
after the summation over A^. . In the first case we find zero for odd values of A^ while for 



r N 

/ VVe-x.Y>{—i:T{mV-^e-KY>{-id/Nf)+mVe-x.Y>{ie/N})), (3.10) 
J dot v=i 2 A 



even values of only the term Nj^ = = N/2 contributes in which case the a integral 

becomes soft. However, in the second case the interference between contributions of all 

different topologies results in the integrand (2cosa)^ which allows us to perform the 

integral by a saddle point approximation in order to obtain the leading order contribution 

in The saddle points are located at a = and a = tt. It is at this point that the 

distribution function fi{N) plays a role: the contribution of the latter saddle point, (—1)^, 

can be ignored after the summation over A^, whereas the contribution at a = yields 

a overall factor jJ,{N) that does not contribute to the m dependent part of the partition 

function to be discussed below. 

At the saddle points in A and a, the [/—dependence can be absorbed into V. The the 

[/—integration yields a finite irrelevant constant. Since the phase expia has already been 

extracted the remaining integral over V is over SU{Nf) instead of U{Nf). We treat m 

as a small parameter and expand the determinants up to first order in m. At the saddle 

point a = the result for m dependent part of the partition function is 

Z{m,0) r .N. 
Z{m = 0,9 

which coincides with the result derived by Smilga and Leutwyler [0] for the QCD partition 
function using chiral perturbation theory. 

The value of S at 6' = can be obtained from eq. (2.4) 

j:(e = 0) = hm hm < Tr(V + V~') >zie=o) • (3.11) 

m.~>0 N~^oo ZX J\ f 

For Nf = 1 the integral over V is absent, and the sign of the quark condensate is inde- 
pendent of the sign of m. For more than one flavor the order of the limits allows us to 
perform the V integral by a saddle point approximation. In the case of equal positive 
masses the saddle point for 6' = is at V = 1 which allows us to identify the parameter 
A as 

' = W^y 

which completes the calculation of the partition function. 

As observed in ||^, for more than one flavor the value of the condensate depends on 
the sign of the quark mass. For example, in the case of two flavors with equal negative 
masses and 6 = 0, the saddle point is at \^ = —1. 
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It is also possible to introduce a complex mass in eq. (2.2) with the mass in the lower 
block equal to the complex conjugate of the mass in the upper block of the overlap matrix. 
In this case the final result for the partition function depends only on the combination 
me*^ and its complex conjugate which makes it possible to derive relations between 
the m and the 9 derivatives of the partition function. 



4. Sum rules for one and two flavors 

Sum rules for moments of the inverse squares of the eigenvalues of the Dirac operator 
can be derived for an arbitrary number of flavors and for any value of the total topological 
charge P]. In this section we only present a derivation for the simplest two cases of one 
flavor and of two flavors with equal masses. Proofs of the general results can be found in 

i- 

For Nf = 1 there is no integration and the result for the ratio of the massive and 
massless partition function is particularly simple 

ZNf=i = exp(A^mS COS 6^), (4.1) 

which was first obtained in 0. 

In the case of two flavors the integral over SU{2) can be performed for an arbitrary 
mass matrix. Here, we restrict ourselves to the simpler case of equal quark masses. In 
this case the ratio of the massive and the massless partition function reduces to 

Z]\f=2= / — sm — exp^iA/mL cos — cos -), (4.2) 
JO TT 2 2 2 

where the factor sin^ | results from the invariant measure of SU{2). The integral is 
elementary and results in 

H2,nN^cosJ^^ ,4.3) 

mNTi cos 2 

The partition function can also be evaluated by writing the fermion determinant as a 
product over eigenvalues. For Nf = 1, we obtain sum rule 

2 

exp{EmN cosO) =< m''ll'{l + —) >z{m=o), (4.4) 

n 



where the average is with respect to the partition function Z with m = 0, which also 
includes the weight factor H'-^n involving the product of the nonzero eigenvalues of the 
Dirac operator. The exclusion of zero eigenvalues in the product is denoted by a prime. 
For two flavors a similar sum rule can be derived 

n (1 + 4m^^^^cos^ I ) ^ ^^^^^^^ (4 5) 

Here, the X\n is the n'th zero of the Bessel function Ji. By expanding both the r.h.s. 
and the l.h.s. of eqs. (4.4) and (4.5) in powers of m we obtain sum rules for the inverse 
moments of the eigenvalues A^. Additional sum rules for N f = 2 can be derived from 
the expansion in powers of both m„ and rrid of the general expression for the partition 
function. 

The sum rules (4.4) and (4.5) are for a fixed value of 6. However, we will investigate 
the validity of the Leutwyler-Smilga sum rules in an interacting instanton vacuum where, 
for technical reasons, the total topological is taken to be zero. The corresponding partition 
function Zq is given by the Fourier transform of Z{6), 

Z, = -Z(B). (4.6) 

which can be derived from the the decomposition 

Z = Y^e'-^Z,. (4.7) 

For one and two flavors, the integrals over Q are elementary. The results 

Zo = /o(ArmS) for A^^ = 1, (4.8) 
Zo = lliNmT) - ll{NmY.) for = 2, (4.9) 

were first derived by Leutwyler and Smilga 0. Sum rules are obtained by equating (4.8) 
and (4.9) to the r.h.s. of (4.4) and (4.5) restricted to zero topological charge, respectively. 
For Nf = 1 one finds 

2 

771 

loiNmE) =< n(l + -tt) >.=o, (4.10) 

n 

and for Nf = 2 we have 

2 

777 ^ 

P.iNm^) - /^(iVmS) =< n(l + >.=o ■ (4.11) 
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The expansion of both sides of this equation in powers of mN yields sum rules for moments 
of the inverse eigenvalues of the Dirac operator. The first two sum rules are 

<E]^>-o = J^, (4.12) 

1 1 1 

(4.13) 

As was shown by Leutwyler and Smilga, both sum rules hold for arbitrary Nf and also 
be worked out for arbitrary topological charge. The only modification is to replace Nf 
by Nf + In the case of two or more fiavors additional sum rules can be derived by 
varying the quarks masses independently 0. We only quote the result 

1 _ 

Again, the result for arbitrary topological charge is obtained by replacing Nf by Nf + 
Formally, this sum rule diverges for Nf = 1. That this is indeed the case can be understood 
from the behavior of the spectral density at small A. For = we expect p(A) ~ 
A^/^^, where the factor A^-'' originates from the fermion determinant and the factor A^ 
from the Jacobian of the unitary transformation that diagonalizes the overlap matrix. 
The above sum rules diverge for Nf = 0. Consequently, we expect that in this case, and 
therefore in the quenched approximation, that the thermodynamic limit of the spectral 
density near zero is approached in a completely different way. 

The sum rules (4.12-14) can be expressed in the microscopic spectral density by writing 
the sum over the eigenvalues as the integral 

V-Ap = !—ps{x), . (4.15) 

The second sum rule expresses an integral over the correlation function < ps{x)ps{x') >u=o 
in terms of integrals over the average microscopic spectral density < ps{x) >u=o- 

Finally, let us write down a sum rule for nonzero values of mN. In the case of one 
fiavor one obtains from the infinite product expansion for the Bessel function Iq, 

^ mNTi^ r°° mN 

\xl^^m^N'^Y? ^ Jo '^"^^^^^'^^l-.ivflxed >-=o 3-2 + ^2^2' (4.16) 
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where xon is the n'th zero of the Bessel function Jq. One is tempted to invert this equahty 
in order to obtain an analytical answer for the microscopic spectral density. Since pg 
depends on mN there is no unique solution, but a trivial solution can be written down 
readily, 

P5(x) = E^(^-^)- (4-17) 

Because there is no reason to believe that the dispersion (4.13) is zero, we do not expect 
that this is the exact microscopic level density. 

In appendix A it is shown that for one flavor sum rules can also be derived at finite 
A^. These results provide an explicit proof that the thermodynamic limit has been taken 
correctly. 



5. Sum rules in the instanton liquid 

The partition function (2.1) can be approximated semiclassically by a liquid of in- 
stantons. Instead of averaging over all gauge field configurations, we average over the 
collective coordinates of the instantons only, whereas 1-loop quantum fluctuations about 
the instantons are included in the measure. The action in (2.1) is the instanton action, 
which also includes the interaction between instantons. We use the so called streamline 



20| , pi| interaction supplemented by a core in order to stabilize the instanton liquid. The 



fermion determinant is calculated in the space spanned by the fermionic zero modes with 



overlap matrix elements that can be derived from the streamline configuration [22|. More 



details on the above instanton liquid model can be found in |T3[ . 

The numerical simulations were carried out for a liquid of 64 instantons in a Euclidean 
space time volume of (2.378)^ x 4.756 in units of A~^. Averages were obtained from 
125 statistically independent configurations. Our main results are presented in Table 1. 
Calculations were done for one, two and three massless flavors (see heading). The first 
row shows the smallest eigenvalue and its average level motion (between brackets). The 
condensate S (row 2) is obtained from an extrapolation of the spectral density to A = 0. 
The results for the sum rule (4.12) (row 3) are compared to the analytical result (row 4). 
We find complete agreement inside the error bars. The remaining rows involve numerical 
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and theoretical values for the quantities Si and 5*2 defined by 

< (E p-)^ >i/=o - < E 3^ >j/=o ^ <(Ep-)^>!^=o , , 

Q n ^^n Q ^^n l K \ 

•^1 - ^ 1 ^2 ' - ,^1.2 • ^^■^) 

Also in this case we find excellent agreement with the theoretical values 



Qth _ ^^f qth _ / /cr n\ 



where the derivation for 5*2 does not hold for Nj = 1. We consider this a strong argument 
in favor of the universality of the sum rules. 



6. Discussion and conclusions 

The main conclusion of this work is that the low energy chiral limit of the QCD 
partition function as derived by Leutwyler and Smilga coincides with a large limit of a 
random matrix theory. Therefore all conclusions of [0 also pertain to the chiral random 
matrix theory proposed in section 2. The random matrix theory resembles the partition 
function of a gas of instantons. The main difference is that instead of inheriting its disorder 
from the distribution of the collective coordinates, the overlap matrix elements are taken 
to be independently distributed according to a gaussian random matrix ensemble. 



From random matrix theory it is know that |23| microscopic level correlations such as 
for example the nearest neighbor level spacing distribution or the variance of the number 
of levels as a function of the average number of levels in a given interval are universal. 
They only depend on the symmetries of the Hamiltonian but not on the details of the 
matrix elements. The spectral sum rules put forward by Smilga and Leutwyler involve the 
microscopic level density. Since they can be obtained from a chiral random matrix theory 
we expect that the sum rules are also universal, and that they do not depend on the details 
of the low-energy structure of the theory but only on its symmetries. In particular, they 
depend on the number of fiavors and, for example, diverge in the quenched approximation. 
This conclusion has been confirmed by simulations of an instanton liquid with overlap 
matrix elements that differ strongly from those of a gaussian random matrix ensemble. 
Inside the statistical uncertainty of the calculations all sum rules were reproduced. The 
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physical interpretation is that the way the thermodynamic hmit of the spectral density 
near zero is approached is universal. However, its value, i.e. the chiral condensate, depends 
on the details of the theory. 

The advantage of a random matrix model is that it is amenable to powerful mathemat- 
ical techniques [^, Q that makes it possible to obtain explicit analytical results for e.g. 
the microscopic spectral density, which, by virtue of the universality arguments discussed 
above, also hold for the QCD partition function. Work in this direction is under way. 
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Appendix A 

In this appendix we derive the sum rules (4.12) and (4.13) for one flavor without 
relying on a saddle point approximation. For A^_|_ = = N/2 and Nf = 1 the partition 
function simplifies to 

/■__JV_JV N 

Z = d(jda{a + ia — m) ^ {a — icr — m) ^ exp(— — — -(o" + ia){(T — icr)), (A.l) 
J IX 

where an irrelevant overall constant has been suppressed. At finite A the pre-exponential 
factors can be expanded as a binomial series which provides us with an expansion in 
powers of m. The coefficients are elementary integrals, and for the m— dependent part of 
the partition function we find 

^--^ = 1 + + 1 +■■■• A.2 

Z(0) 4 64 ^ A^ ^ ' 

It should be noted that no approximations have been made. On the other hand, the 
fermion determinant can be written as a product over the eigenvalues which leads to the 
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expansion 

where the average is with respect to the massless partition function. By equating the 
coefficients of the powers of we obtain sum rules for the inverse powers of the eigen- 
values that are vahd for any value of A^. The sum rules (4.12) and (4.13) for Nf = 1 
are reproduced by keeping only the leading order terms in 1/A^. We observe that the 
first sum rule is valid for any value of N. The second sum rule is modified by the factor 
(1 — 2/N). For N — 2we find zero which is correct because in this case there are no terms 
that contribute to the sum n' ^ n in (A. 3). 
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Nf = 3 



Ai 

E 



0.0255(98) 



0.0531(164) 



0.1078(335) 




1.90(5) 
0.89(3) 
0.90(3) 
0.53(3) 
1/2 



1.33(5) 
0.21(1) 
0.22(1) 
0.70(5) 

2/3 



0.90(10) 
0.066(2) 
0.068(7) 
0.83(7) 

3/4 



2.0(4) 



1.20(12) 
4/3 



1.24(12) 
9/8 



oo 



Table 1. Numerical results for the spectral properties of the Dirac operator in the gauge 
field of a hquid of instantons. Prom the comparison of the 3rd and 4th, the 5th and 6th 

and the 7th and 8th rows we conclude that the Leutwyler-Smilga sum rules are observed 
by this model. For the definition and further discussion of the observable, we refer to 
section 5. 
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